HYPOCOERCIVITY FOR LINEAR KINETIC EQUATIONS 
CONSERVING MASS 

-^ , JEAN DOLBEAULT, CLEMENT MOUHOT, AND CHRISTIAN SCHMEISER 

^,^' Abstract. We develop a new method for proving hypocoercivity for a large 

J^ ' class of linear kinetic equations with only one conservation law. Local mass 

conservation is assumed at the level of the collision kernel, while transport in- 
volves a confining potential, so that the solution relaxes towards a unique equi- 
librium state. Our goal is to evaluate in an appropriately weighted L norm 
the exponential rate of convergence to the equilibrium. The method covers 
various models, ranging from diffusive kinetic equations like Vlasov-Fokker- 
Planck equations, to scattering models like the linear Boltzmann equation or 
P^ ' models with time relaxation collision kernels corresponding to polytropic Gibbs 

^^ ' equilibria, including the case of the linear Boltzmann model. In this last case 

^J and in the case of Vlasov-Fokker-Planck equations, any linear or superlinear 

r^ , growth of the potential is allowed. 

-)— > 



> 



1. Method, result and consequences 

1.1. Linear kinetic equations and hypocoercivity. We consider linear kinetic 
ly-N ' equations which can be written as 

^; (1.1) ft/ + T/=L/, f = f{t,x,v), {t,x,v)eR+xR''xR'', 

and describe the evolution of a distribution function f. The transport operator 

C^ I has characteristics given on the phase space M.'^ x M'* by the flow of the Hamiltonian 

{x, v) ^ E{x, v) := - \v\^ + V[x) . 

K^ \ The external potential V — V{x) is a measurable function on M''. The collision 

; I ■ operator L is independent of time t and acts as a multiplicator in the position 

Cu I variable x. The variable v is the velocity. 

We shall consider steady states which are in the intersection of the null spaces 
of T and L simultaneously We shall assume that there exists a nonnegative energy 
profile function F such that, for each fixed value of x, the nuUspace Af{\-) of L is 
spanned by F{x, v) :— T{E{x, v)), so that 

N{L) = {f{x, v) : 3 (f>{x) such that /(x, v) = (f>{x) F{x, v)} . 

Functions in A/'(L) are local equilibria; they depend on x and t. The function F 
is a global equilibrium or global Gibbs state. It is independent of t (stationary) 
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and isotropic with respect to v. Consistently, we shall further assume that L has 
rotational symmetry in v, i.e. Ri, L = L R^, for any rotation operator Ry acting on 
the velocity space. Under the assumption that the support of F is connected, the 
intersection of the null spaces of L an T is generated by F. Assume that F is 
integrable and normalized by 

F dvdx = 1 . 

d 

We shall refer to this assumption as Assumption (HO) and assume that it holds 
throughout the paper, although we shall not specify it explicitly when it is not useful 
for the understanding of our arguments. Under such a normalization condition, we 
shall prove that F is the unique stationary distribution function. Integrability with 
respect to v is an assumption on F, whereas integrability with respect to x requires 
a F-dependent growth of the external potential V. Such a property is a confinement 
condition. 

The one-dimensionality of Af{\-) for fixed x suggests the existence of one local 
(in x) conservation law. We shall therefore assume the local conservation of mass, 
that is 

Lfdv = 0. 
Global mass conservation for solutions of (|l.ip follows: 

^ [[ fdvdx= [[ {L~T)fdvdx = 0. 
For an integrable initial datum 

/(i-0, •,•)-//, 

let M := JJg^^^^j fj dvdx, so that MF is the unique global Gibbs state with 
mass M . In this paper we investigate the asymptotic behavior of the semigroup 
generated by L — T. Our goal is to quantify its stability or, to be precise, to 
determine the rate of convergence of / towards MF as t — >■ cx). Since the equation 
is linear, there is no restriction to study fluctuations around a global equilibrium, 
that is solutions / of p.ip which satisfy 

(1.2) M= ff f{t,x,v)dvdx^ [[ fi{x,v)dvdx^O. 

Notice that distribution functions are usually nonnegative, but fluctuations around 
an equilibrium have to change sign. 

Local mass conservation for / and F imply the identity 

for any function (j) ~ (j){x), thus showing that the left hand side is, at least formally, 
quadratic in the distance between / and the kernel of L. This suggests to introduce 
the space L^{d^) where the measure d/i is defined on the phase space by 

fii) ri IT 

d^i^d^iix^v):^ —- r, (x, w) e M'' X M''. 

b [x, v) 
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We shall denote by (•,•) the corresponding scalar product and by || • || the asso- 
ciated norm, so that (/,.g) = ^j^dy^df gd^i and ||/|p = {f,f). The orthogonal 
projection 11 on the set of local equilibria is denoted by 



n/:=^F, withp/:- / fdv. 

PF jRd 



We also assume that the collision operator is dissipative in the sense that an 'H- 
theorem' holds, i.e. (L/, /) < 0. Since the transport operator T is skew symmetric 
with respect to (•, •), this implies the entropy inequality 

^|ll/irHL./,/)<o. 

Under the normalization condition (11.21) . if the entropy dissipation —{Lf,f) was 
coercive with respect to the norm || • ||, exponential decay to zero as i — ?► cx) would 
follow. However such a coercivity property cannot hold since L vanishes on the set 
of local equilibria. Instead we shall assume that microscopic coercivity holds, i.e. 
there exists a positive constant Am, such that 

-(L/, /) > A„, 11(1 - n) /||2 for all / G L^dfi) . 

The key tool of our method is a modified entropy functional H[/], whose square 
root is a norm equivalent to 1 1 • 1 1 , such that 

|h[/]<-AH[/], 

for an explicitly computable positive constant A. As a consequence, we find an 
estimate of the exponential decay rate of the semigroup. Following the vocabulary 
used in [33l [20l [26] , such a strategy will be called hypocoercivity. 

In some cases, the existence of a spectral gap can be obtained by non-constructive 
compactness methods, see for instance [35] in the case of the linearized Boltzmann 
equation on the torus. For a non-positive closed operator U with a spectral gap 
A > 0, it is well-known, see [27], that there exists a norm equivalent to the ambiant 
norm, for which the semigroup of U -|- A is contractive. However this method is not 
constructive regarding the norm of contractivity and gives no estimate on A. In 
our approach, under assumptions specifically adapted to kinetic theory, we are able 
to construct an explicit Hilbert norm which is equivalent to the standard norm of 
L^{dfi) and to estimate A. 

Various results related to hypocoercivity have recently appeared, on large time 
estimates: [13l O [8] [9]; based on hypoellipticity: [211 Ull 122]; on hypocoerciv- 
ity itself: [33l (20] [26] ; on applications of the so-called kinetic-fluid decomposition: 
[T^ [HI [TBI [T71 [T51 ^, '3D', '31] ; on hyperbolic estimates based on micro- macro de- 
compositions: [23l[24, 25, 34:. Some of the results of this paper, namely Theorems[9] 
and ll4[ have been announced in [11] without complete proofs. 

Our purpose is to establish, in a simplified framework, sufficient conditions for 
proving hypocoercivity for a large class of linear kinetic models confined by an 
external potential, without assuming regularity on the initial datum and valid for 
hypoelliptic kinetic Fokker-Planck equations as well as singularity preserving col- 
lisional kinetic equations. This is the main difference with hypoelliptic methods. 
The method also makes use of a micro-macro decomposition. Accordingly we shall 
split our assumptions into two main requirements: microscopic coercivity as intro- 
duced above, and a macroscopic coercivity assumption, which is a spectral gap-like 
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inequality for the operator obtained when taking an appropriate macroscopic dif- 
fusion hmit that we shall now describe. 

1.2. Formal macroscopic limit. As a motivation for the macroscopic coercivity 
assumption, we recall, at a formal level, the macroscopic diffusion limit procedure, 
which can be seen as intermediate asymptotics governing the long time behaviour 
of solutions. On a large time scale, it heuristically models how local equilibria 
relax towards the global Gibbs state. Since the macroscopic flux of the equilibrium 
distribution vanishes, i.e. J^g^ v F dv = 0, the appropriate macroscopic rescaling of 
the solution of (|1.1|) is given by 

f^{t,x,v) = f{e'^t,ex,v) , 0<e<l, 

which is known as the parabolic rescaling. Assuming that the potential V is rescaled 
accordingly, we obtain the singular limit problem 

as e — >■ 0. The assumption lime^o f^ — /" leads to L/° = and, thus, /° = 11/°. 
The identities 111 = Ln = nTn = imply the relations 

edtf +Tf ^\-R' , dtUf +UTR' = 0, 

with i?^ := i (1 — n) /^. Assuming formally that lime^o ^^ — R'^, the first equation 
can be solved for e = with respect to i?°, giving 

i?° = JT/° = JTn/°, 

where J denotes the inverse of the restriction of L to the orthogonal complement of 
its null space. Note that the inhomogeneity T 11 /"^ satisfies the solvability condition 
XITn/" — 0. The second equation becomes 

(1.3) ^tUf = i^ur^i^u)f, 

where the superscript * denotes the adjoint operator with respect to (•,•), and the 
skew symmetry of T has been used. A straightforward computation shows that this 
is equivalent to a drift-diffusion equation for the macroscopic density p'^ = pfo : 

(1.4) atp° = V, • [p^aV,(|^)] =V, • [v,(f7pO)+7pOV,F 

Here a is scalar due to the rotational symmetry of L, 

If 1 

PfO' — —-: v-]{vF)dv and ^'^xV— "^xiPFC^)- 

d J^d pF 

The operator J being negative definite on (1 — H) L'^ldp.), a{x) > for all x E M'*. 
In the two following important cases, the macroscopic transport coefficients 7 and a 
have particularly simple expressions. 

Case (CI). When r(s) = e~*, the global Gibbs state is a Maxwellian, or Gaussian 
function, which factorizes as 

F{x, v) = pf{x) M{v) , with pp = j _y and VJl{v) = ^.^ . 

Jjjd e ax (Ztt) 

Notice that the separation of position and velocity variables is a characteristic 
property of Maxwellian functions. Both coefficients 7 and a are constant, equal to 
2 /jjd V ■ J{v Wl) dv and p° solves the Fokker-Planck equation 

dtp° ^aVx-{Vxp'' + p°VxV) . 
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Case (C2). The collision operator L is, for fixed x, a time-relaxation operator 
onto span{i^}, i.e., 

L = n-i. 

In this case J = —Id holds, so that 

pp(T = mp:=— / \v\'^ F dv , 

and, since Vx{pf o') ^ ^ Pf ^xV because 

i/ \v\^r'{E{x,v))dv = - f vVyFdv^-pF, 

the macroscopic limit equation reads 

9tP° = V, •(V,(ap°)+pOV,l/) . 

The intersection of both cases, (CI) and (C2) i.e., L = 11 — 1 with li f = pf 9Jl, 
gives CT = 7 = 1. This is the linear BGK case, which has been considered in [llj . 
In both cases, (|1.3p can be rewritten as 

9tn/° = -ao(Tn)*(Tn)/o 

for some positive constant <tq, with a^^ a in Case (CI) and (Jq = 1 m Case (C2). 
With Assumption (|1.2p on the initial data, we expect decay to zero of the solution. 
Under a macroscopic coercivity assumption, namely (H2) (see below), which is 
equivalent to a Poincare inequality (see Lemma lTTS)) .the decay of IT /° is exponential. 

1.3. Method and main result in an abstract setting. We start with the basic 
assumption that L and T are closed linear operators on an Hilbert space H, such that 
L — T generates the strongly continuous semigroup e'^'-^^^* on H. The orthogonal 
projection on the null space Af{\-) of L is denoted by 11 and 'D{L) is the domain 
of L. We assume that the restriction of L to Af{L)-^ is coercive. More precisely, our 
first assumption is: 

Assumption (HI) (microscopic coercivity): The operator L is symmetric and 
there exists Am > such that 

- (L /, /) > A,„ II (/ - n) /f for all f eV{L). 

Motivated by the results of Section 11.21 coercivity of the transport operator is 
required, when acting on A^(L): 

Assumption (H2) (macroscopic coercivity): The operator T is skew symmetric 
and there exists Aj\/ > such that 

||Tn/f > AmIIH/IP for all f en such thatUf eV{T). 

Inspired by j20| . we introduce the modified entropy 

H[f] := i II/IP + e (A/,/) , with A := (l + (Tn)*(Tn))-\Tn)* . 

The constant e > will be chosen below. A straightforward computation for a 
solution / of ()1.1|) . now considered as an abstract ODE, gives 

|h[/] = -D[/] 
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where the dissipation of entropy functional is given by 

D[/]:--(L/,/)+e(ATn/,/)+£(AT(l-n)/,/)-£(TA/,/)-£(AL/,/). 

By (HI), (H2), and by ATn = (1 + (Tn)*(Tn))-i(Tn)*(Tn), the sum of the 
first two terms in D[/] is coercive: 

-(L/,/) +e(ATn/,/) > min{A,„,^i^} Wf^ 

For the completion of our program, we need to show that H[/] is equivalent to ||/||^ 
and to control the last three terms of D[/]. Part of this can be carried out at the 
abstract level under the following additional assumption: 

Assumption (H3): 

nTn = o. 

Lemma 1. Let Assumptions (HI)— (H3) hold. Then the operators A and T A are 
bounded, and for all f ^ %, 

(1-5) ||A/||<i||(l-n)/|l and |1T A/|l < |1(1 - H) /|1 . 

Proof. The equation A f = g is equivalent to 

(Tn)7 = 5 + (Tn)*(Tn)5. 

Writing this as g = UT^Ug -UT f proves A = HA and, thus, TA/ = TUg. 
Taking the scalar product of the above equality with g and using (H3), we get 

llgf + ||Tn.g|p = (/,Tn.g) = ((l-n)/,Tng) 

< ||(i-n)/||||Tn5||<i||(i-n)/|p + ||Tn.9|p, 

which completes the proof. D 

The boundedness of the remaining terms in D[/] has to be proven case by case. 
We shall therefore assume it in the abstract setting. 

Assumption (H4) (Boundedness of auxiliary operators): The operators AT (1 — 
n) and A L are bounded, and there exists a constant Cm > such that, for all f Cz H, 

||AT(i-n)/|| + ||AL/||<CM||(i-n)/||. 

Theorem 2. Let Assumptions (HI)— (H4) hold. Then there exist positive con- 
stants A and C , which are explicitly computable in terms of Xm, Xm , and Cm, such 
that, for any initial datum fi^T-L, 

||e*(L-T)//|| <Ce-^*||M|, Vi>0. 

Proof. The first inequality in (jl.Sp implies 

(1-6) \{l-e)\\ff<y\[f]<\{l + e)\\f\\\ 

For any e £ (0, 1), H[/] is equivalent to ||/|p. The second inequality in (|1.5p and 
(H1)-(H4) imply 

D[/] > A„j|(i-n)/|p + ^|in/f-e(i + CM)ll(i-n)/||||/|| 

> [A„, - e (1 + Cm)(1 + M 11(1 - n) ,/lP + e [^ - (1 + Cm)I\ ||n/f 
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for an arbitrary positive S. By choosing first S and then e small enough, a positive 
constant k can be found, such that D[/] > k ||/||^. Using (|1.6[) . this implies 

|h[/]<--^H[/], 

at 1 + e 

for / = e'^""^^^//, completing the proof with A = k/(1 + e) and C = \/ jz^- □ 

Let us conclude this abstract approach by some comments. First of all, our proof 
is constructive: H is an explicit Lyapunov functional and A can be computed. The 
work of F. Herau has been a crucial source of inspiration for our method. In "^ , he 
deals with the linear time relaxation collision kernel corresponding to Maxwellian 
Gibbs states, in case of a confining potential V growing at most quadratically 
at infinity and such that the associated Witten Laplacian satisfies a spectral gap 
inequality. In our approach, we are able to relax some of these assumptions. See 
Theorem P 

Our results apply to various Fokker-Planck and Boltzmann models. We shall 
compare applications of Theorem [2] to previous results in Section 3. Only [2 6) and 
[33] deal with abstract results like the ones of Theorem [2] Ours are more general 
than the ones of [5S] since we deal with a general confining potential. In [5B], the 
problem is indeed set on a torus, a setting to which our method can be adapted 
without any difficulty. It is also more general than in [33_ since we deal not only 
with Fokker-Planck type operators, or operators in Hormander form in the words 
of [33], but also with non-local integral collision operators, like in [26]. Last but 
not least, our results are also stronger than those in j^S] and [33] in the sense 
that we construct a zeroth order norm of hypocoercivity, which is equivalent to 
L^ and not iJ*^, for some fc > 1. However, our results are weaker than those in 
[25] at least in one aspect: we only deal with models with 1-dimensional space of 
collision invariants, whereas, in |26j . any finite dimension is allowed. In principle, 
our approach can be extended to such a situation, which is the purpose of a current 
research project [12] . 

1.4. Hypocoercivity for a toy problem. To illustrate the fact that our formal 
setting applies to other models than the kinetic equations of Section 11.11 we intro- 
duce the following toy model, which captures very well the essential features of our 
hypocoercive approach. We consider a one-dimensional Cattaneo model introduced 
in [6], which can be written as a kinetic model with only two velocities w = ±1, and 
where L describes a switching process between the two velocities without preference 
for one of them. As a further simplification we replace the confining potential by 
a periodicity assumption, where x varies in a one-dimensional torus. The model 
equations are 

for the distributions f^{t,x) of right- and left-moving particles, periodic in x with 
period 2 tt. 

The interest of such a model is that it gives an application of our hypocoercivity 
method in a discrete setting, or even for a finite dimensional ODE version of it, 
if we truncate the Fourier sum in the x variable and keep only a finite number of 
terms. 
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Initial value problems can be solved explieitly by Fourier decomposition. Intro- 
ducing the total density p = f^ + f^ , the total flux j = f^ — f^, and their Fourier 
representations 

pit,x) - ^Pfe(t) e^"^ , j{t,x) - 5^ jfe(t) e^"^ , 
fcgz fcez 

leads to real ODE systems for [/,. = {^D := (f^g^;J) and Uk :- CZul))' 

(1.7) ^+T,Uk = LUk, 

at 

where the skew symmetric matrix T^. := (^ ^'^) represents the transport operator, 
L := (o _?i) represents the collision operator acting only on the microscopic com- 
ponent j, and Uk solves an analogous system with T^ replaced by — T^. Eq. (|1.7p 
is linear, and it is elementary to check that the eigenvalues of L — T are given by 
^o,± '■— 0,-1 and Afc^± := (—1 ± i^Ak'^ — l)/2 if fc 7^ 0. All solutions converge 
to an eigenstate of the zero eigenvalue: C/q = (po,0) and Uk = ioi k ^ 0. The 
convergence is exponential with its speed determined by the spectral gap 1/2. 
For fc 7^ 0, we can compute the entropy dissipation as 

so that it is clear that no exponential decay directly follows, since the right hand 
side is not coercive and there is an unbounded increasing sequence (in)neN such 
that Vk{tn) = 0. Note that microscopic coercivity holds with A™ — 1- 

With n = (1 0) and J = -Id, we find that (Tn)*J (TH) = ( "fc" 0), thus giving 
for the macroscopic diffusion limit du^/dt = —k'^u^, and showing also that macro- 
scopic coercivity holds with Am = 1- According to the strategy of the Section [T73l 
for fc 7^ we introduce the modified entropies 

Hfc W = \ \Uk{t)? + e Y^^ uk(t) vk(t), t>0. 
Observing that, for fc > 1, 

'2{^-^TTw)\Uk\'<^,{l-ej^)\Uk\' + ^,j^\uk+Vk\' = Hk, 

using supj.>]^ ,',i = ^ < I7 and performing a similar computation for k < —1, we 
finally get 

Ul-e)\Uk\^<Hk<'^{l+e)\Uk\'. 

Hence, for any e G (0, 1), Hfc(t) decays exponentially if and only if \Uk{t)\'^ decays 
exponentially as well. Obviously, we have 

1 k^ 
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which makes it easy to compare Hfc with ^ H^, given by 

for any A G (0, 1). If £ e (O, ^\i,i ), the coercivity constant 

K := min {§ (1 — A^), 1 — e — g^} 
is positive and 

|H.<-.|t/.P<-^H.. 

This impHes |(7fc(i)| decays hke e~''*/^^+'^^. We may observe that 

K min{f,l — e} 1 
< ^-^ < 



1 + e 1 + e - 5 ' 

thus showing that the method is not optimal, in the sense that it does not give the 
exact decay rate, 1/2, even when refining the above estimates and computing k for 
each k. 

1.5. Application to kinetic equations. Let us apply the abstract procedure of 
Section [T751 to the setting of Section [TTTl Thus, we set 

PF 

where the potential V is given as well as the energy profile F. We recall that the 
unique global equilibrium is F{x,v) = T{\v^\/2 + V{x)), x,v £ R"*. For such an 
equilibrium distribution, define the velocity moments up to the fourth order by 



PF 



I Fdv, mF-=^ [ \v\^Fdv, Mf:= [ \v\^F dv 



and assume that they are measurable functions of a;. We consider the Hilbert space 
Ti, — {f E L^{d^) : JJ^ay^jg^d f dvdx — 0}, with d^{x,v) = dxdv/F{x,v). The 
collision operator L remains unspecified at the moment, so that we shall defer the 
discussion of the microscopic coercivity for a while. A simple computation with 
u = Pf/pF shows that the macroscopic coercivity assumption is equivalent to a 
weighted Poincare inequality: 

Lemma 3. Assumption (H2) holds if and only if 

(1.8) / IV^:^! niF dx > X]\f / u^ppdx 

for any u G L^ipF dx) with 'S/xU G L^lmp dx) such that L^ u pF dx — 0. 

In case of kinetic equations. Assumption (H3) is a consequence of the compu- 
tation 

JUf ^FvVxUf, 

with Uf :— pf/pF, and of the observation that the right hand side is an odd 
function of w, whose mean value is zero. In other words: The macroscopic flux of 
the equilibrium distributions vanishes. 
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Concerning the Assumption (H4) , we remark that boundedness of A L is pos- 
sible even for unbounded colhsion operators L (see Section [2]) . The boundedness 
assumption on A T (1 — H) can be interpreted as an ehiptic regularity result for: 

(1.9) U Va; • {mp Vj;U) = w . 

PF 

Lemma 4. // there exists a positive constant C such that 

(1-10) \\^lu\\L'^(Mpdx) < C\\w\\L2(^p^dx) 

for any w G L^{pf dx) and for any solution u G L^{pf dx) with "S/ xU £ L'^imp dx), 
then the operator AT (1 — 11) is bounded on %. 

Proof. The operator AT (1 — 11) is bounded if and only if its adjoint 

[AT(i-n)]* = -(i-n)T2n[i + (Tn)*(Tn)]-i 

is bounded. If .g = [1 + (Tn)*(Tn)]- V, then 

[AT(l-n)]*/ = -(l-n)T2ng with .g + (Tn)*(Tn)g = /, 

where the latter implies (jl.9p for u — Ug ^ Pgl Pf a-nd w = Uf = pf/pp. Then 

T^ng = F V ■ V x{v ■ V xu) + FV xV ■ V xu 

results in 

[AT (1 - n)]*/ = -(1 -IV)[Fv ^x{v ■ Vxu)] . 

This implies that for some positive constant C, we have 

||[AT(1 - n)]*/ll <\\Fv V,(« • V,u)|| <c\\Vlu\\L2iM,dx) , 
which completes the proof using ||u/||l2(pp ^j.) = ||n/||. D 

2. A FRAMEWORK FOR THE ELLIPTIC REGULARITY ESTIMATE 

Our goal is to give conditions on V which are sufficient to establish the existence 
of a positive constant C as in Lemma |4l In the applications considered below, the 
combination of weights Mp pp/m\ is constant. This motivates the notations 

wl'.^pp, wf -.^^ {jfYwQ with i = l,2, ||u||i := ||uw.,||i,2(Rd) . 

With p = upp, the Poincare inequality in (H3) can then be rewritten as 

(2.1) \\^xu\\1>Xm\\u\\1 

under the zero average condition J^^ u pp dx — 0, and the desired estimate (|1.10p is 

(2.2) \\Vluh<C\\uf\\o 
for the solution of 

(2.3) wlu~Vx- {w\Vxu) ^wluf . 

Roughly speaking we just have to prove (L^ -^ iJ^)-regularization for a second 
order elliptic equation. However, different norms have to be taken into account. 
The result can only be shown under certain assumptions on the weights, which will 
later be translated into assumptions on the confining potential: 

(2.4) 3ci > , C2 e [0, 1) such that - wl A2;(logwi) < ci Wq + C2 |Va;Wi|^ , 
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(2.5) 
3 C3 > such that — \W^W\<C3 (l + L^!^ 



(2.6) 



(2.7) 



3 C4 > such that 



IIW^IIo 






< Ci 



with W := \ll 

VajWil 



WxWi 



Wo 



Wo 



W"^ Pf dx < oo . 



Note that a condition on the third weight function W2 could be deduced from (|2.4p - 
(|2.5p - (|2.6p since any two of the weight functions determine the last one. The goal 
of this section is to prove the following i/^-regularity estimate. 

Proposition 5. Let (PTT]) . ([^ . (^3)) . ([^ . and (P77l) hold. Then the solution u 
of dlH) satisfies ((2:21) . 



By Lemma m this shows that the operator AT (1 — 11) is bounded. 

2.1. Improved Poincare inequalities. We start with an improvement of the 
Poincare inequality (|2.ip . 

Lemma 6. Let (12.111 and (12.41) hold. There exists k > such that 



for any u £ L [pp dx) with I u pp dx — 



(2.8) \\VccU\\l>K 



V^wi 



Wo 



Proof. With the identity wiV^u — ^/^{wi u) — uV^-wi, the inequality 
IIVkUJI]^ > u jVajWil dx — 2 uWxWi ■ Wxiuwi) dx 

u^ w\ A3, (log wi) dx 



WxWi 



Wo 

is easily derived. Now (|2.ip and (|2.4I) imply 

VxWi 2 



|V,u||f >(1-C2) 



Wo 



^liv.«ll?. 

A A/ 



This completes the proof with k = Xm (1 — C2)/{Xm + ci). 

Lemma 7. Let ((2T|l . (ITIll . (123]) anrf (l2Jll /loM. T/iere eiiste k' > such that 



D 



\WVxu\\l>n' 



Wu 



VxWi 



Wo 



for any u £ L"^ {pp dx) with / u pp dx — 



Proof. We apply (|2.8p with u replaced by {uW — u) with u := J^^ uW ppdx. We 
recall that J^^ pp dx — JJ^^a^^a F dvdx = 1. We thus obtain 

k||(uH/-?i)h/||^ < \\VxiWu)\\l. 

By expanding the left hand side, we get 

KllVK^ullg < \\\/x{Wu)\\1 + 2k f W^uw^dx ||VK||o||m||o- 
Using 

W'uWppdx<^\\W'u\\l + ^\\W\\l 



12 J. DOLBEAULT, C. MOUHOT, AND C. SCHMEISER 

with a := 2 ||VF||o ||^t||o, we obtain 

2^i [ W^uwldx\\W\\o\\u\\o < na(^\\W^u\\l + ^\\W\\l 

= !l\\W'u\\l + 2K\\Wro\\u\\l 
On the other hand, we can also expand the square in ||Va;(W^u)||^ 



\V,iWu)\\l<2\\WV,u\\l + 2 



uV^W- 



Wo 



and by (P3|) . get 



uVo.W 



Wo 



<ci 



(l + ^)|[<2c^liM^"l 



u 1 + 



Collecting all terms, we finally end up with 

K \\W^ u\\l < 2 \\WV^u\\l +AcI\\Wu\\1 + ^\\W^u\\1 + 2k \\W\\t \\u\\l . 
which, using (|2.1I) . (|2.8p and W > I establishes the inequality: 



3 



^\\W'uro<{2 + A'''^ 



2 — \\W\\t)\\WW^u\\l 
Am ' 



2.2. The regularity estimate. Now we start working on Equation 
standard energy estimate gives 

(2.9) ll"l!o+l|V.ii||?<h/||^ 



D 
The 



With W = ■\/l + iV^jtuiP/wQ , Lemma [5] leads to the improved L^-estimate 

\\uW\\o<-\\uf\\o. 
Lemma [6] can also be used to get an improved _ff ^-estimate. 

Lemma 8. Let (PTT]) . ([^ . (^3)) and (P?7)) hold. Then any solution of (^3)) such 
that Lj u pf dx — Q satisfies 

\\WVM\i<C\\uf\\o. 

Proof. Multiplication of (|2.3I) by u W"^ and integration gives 

(2.10) \\uW\\l+\\WV^u\\l<\\tiW^\\o\\uf\\o~ [ wlu^.,u-y^{W^)dx. 

Since W^ <1 + W \\/xWi\/wo, ((2?T|) and Lemma [7] imply 

(2.11) II" W^'llo < 2 (1/Am + 1/k') ||Vt^ V,u||? . 
The integrand in the last term above can be estimated by 

Wo \u\ wi \y^u\ -^ — — - < 2 V2c3 Wo \u\ wi |V^m| W"^ 

Wo 

with the help of (|2.5I) . so that the integral is bounded by 

2 V2c3 ||uVK||o |1M^V,m||i < 2 V2^ IJM/llo ||W^V,u||i . 

K 

Combining our results gives 

||M^V,u||?<C||u/||o||W^V,u||i 
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with C = ^2(1/Aa/ + 1/k') + 2 \/2 ^, thus completing the proof. D 

Proof of Proposition O We follow the standard procedure for proving _ff ^-regularity 
of the solutions of second order elliptic equations with L^ right hand sides: mul- 
tiply (j2.3p with Vj; • {W xU wf / Wq) , and integrate by parts twice. We also use the 
consequence W2 = wl/wo of the relations between the weights: 



/ Wo {u — Uf) W2 Aa;U dx + Wq {u — Uf) \7xU ■ Vk — 2 I '^^ 

wlVxuVluVxi^] dx- / ^VxuVluV4wl)dx 

X^Q / jRd Wq 



2 



iv^ll^ 



VxU ■ Va; ( ^ ) [VxU ■ Vx (wi)] dx 

= h + l2 
The first integral is easily estimated: 

|/i|<(hl|o + ||H/||o)|lV^^|l2<2|lu;|lo|lV^«|l2 
by (|2.9p . For the second integral we use (|2.6p and Lemma [U 

1/2 1 < 2C4 / Wn (u — Uf)wi\\7xU\ — dx 

Jr-' Wo 

< 2 C4\\u - Uf\\o\\W V xu\\i <4c4C'||u/||o. 
With the third, fourth and fifth integrals we proceed similarly: 

I/3I < 2C4 I Wi\WxU\w2\'^lu\^-^^dx 
Jri Wq 

< 2c4|lW^V,u||i||V2u||2 <2c4C||u/||o||V2u||2, 
I/4I < 2 f wi\Vxu\w2\'^lu\^-^^dx<2c4C\\uf\\o\\Vlu\\2, 

jR'i Wo 

\h\ < 4/ u^l\Wxuf\Wx{^)\^-^^dx 
jRd ' ^Wo" Wo 

<4c4 / wl\Vxuf^-^^^^dx<4ci\\WVxu\\l<4c4C\\uf\\l. 

JR'^ Wo 

The combination of our results gives 

\\Vlu\\l<K\\uf\\o[\\uf\\o + \\\7lu\\2) 
for some explicit constant K > 0, which completes the proof. D 

3. Maxwellian equilibria 

When the local equilibrium is a Maxwellian distribution, the global equilibrium 
has the form 

(3.1) F{x,v)^pF{V{x))9n{v), with pf{V) ^ e-"" and 9JT(z;) = — -^ . 

In this framework. Assumption (HO) is a consequence of 

Assumption (HO.l) The external potential V € C^(R'') is such that e"^ e 
L\dx). 
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As far as the macroscopic coercivity condition (H2) and the boundedness of 
AT(1 — n), i.e. the first part of (H4), are concerned, no further details of the 
coUision operator are required. Consider first the issue of equivalent conditions for 
(H2). 

With w :— pe^/^ the macroscopic coercivity condition (|1.8p is equivalent to 






w dx , 



under the orthogonality condition L^ we^^''^ dx = 0. The first eigenvalue of the 
Schrodinger operator % := — A + i iV^jT^p — ^ AF is zero. It is non-degenerate, and 
the corresponding eigenfunction is w = e^^/^. According to [28], Inequality ()1.8|) 
holds if and only if the lower end of the continuous spectrum of Ti. is positive, that is 

Assumption (H2.1) liminf|^|^oo {\'^ xV\'^ -2AV) > 0. 

As a consequence, macroscopic coercivity holds if IS.V is negligible compared 
to |Va;T^p as \x\ —> cx), and if liminf|2,|^oo iVa;!^! > 0. An example of such a 
potential is V{x) = (1 + jxp)'^ for some /3 > 1/2. See for instance [33l A. 19. 
Some criteria for Poincare inequalities, page. 137] for an elementary proof if 
lim|^l_>oo (I Va;Fp — 2 AV) = cx), and [1] for some recent considerations on Poincare 
inequalities when e~^ is a probability measure. 

Since all three weights pp, itt-f, and Mp are constant multiples of e^^, the 
framework of Section[2]can be used for the boundedness of A T (1 —11) . Assumptions 
([2^ . ([231) are satisfied if 



Assumption (H4.1) There exist constants ci > 0, C2 G [0, 1), and C3 > 0, such 
that 

A.,V < ci + y \VxV\^ , \VlV\ < C3 (1 + IV.V^I) . 

Assumption (j2.6p holds trivially (since wi/wq — const) ^ and (j2.7p can be trans- 
lated to 

/ I'^xVl^e-^dx <oo, 

which follows from (HO.l) and (H4.1) by 

\W^V\^e-^dx = -/ V^V -V^e-^ dx= \ A^V e'^ dx 



< ci / e-^dx + ^ \\/^V\^e-^dx. 
3.1. BGK operator. For the BGK collision operator 

L =n-i, 

the microscopic coercivity condition (HI) is trivially satisfied with A,„ = 1, and, 
since L is bounded (by 1), the boundedness of A L follows from Lemma [TJ 

Theorem 9. Let L = H— 1, and let the external potential satisfy (HO.l), (H2.1), 

and (H4.1). Then solutions of (jl.ll) with initial data in L^{d^) decay exponentially 



to the global equilibrium given by p.l|) . 

This result is an improvement upon the work of Herau !2Q , since the require- 
ments for the external potential are weaker. In particular, Herau's result requires 
potentials with at most quadratic growth at infinity, whereas an arbitrary super- 
linear growth is permitted by (H2.1), (H4.1). 
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3.2. Fokker-Planck operator. For the Fokker-Planck collision operator 

the microscopic coercivity condition (HI) is equivalent to the Poincare inequality 
for the Gaussian measure 93T(u) dv, which satisfies 

2 



(L/,/) - e^ m 



(i-n)/ 



Vt, 2j^ 



> A„,, / eW 971 



(i-n)/ 
m 



dv dx 

2 



dvdx = Xrn\\il~Il)f\\' 



A somewhat surprising fact is the boundedness of A L, although L is an unbounded 
operator. Since A = -(1 + (Tn)*(Tn))-inT and HT/ = -^V^, • j / where jf is 
the flux given by jf = J^a "v f dv, the identity ji / = —jf implies A L = —A, and the 
boundedness of A L is a consequence of Lemma [1] 

Theorem 10. Let Lf — V^ • (Vt,/ + v f) and assume that the external potential 
satisfies (HO.l), (H2.1), and (H4.1). Then solutions of (jl.ip with initial data in 
L^{dfj,) decay exponentially to the global equilibrium given by (|3.1|) . 



The above assumptions are similar to those of [33 , which are weaker than those 
of [21] . Moreover this result is an important improvement compared to |33j as 
involves an L^ setting rather than a H^ setting. Let us emphasize that the latter 
point is not a technical issue and answers an open question raised in |33[ Part II, 
Section 13]). 

3.3. Scattering operators (without detailed balance). Consider a scattering 
operator that can be written as 

(3.2) (L /)(i;) = f [k{v* -^ v)f{v*) - k{v -^ v*)f{v)] dv* , fc > 

where k{v* — > v) denotes the transition probability of changing the velocity v* 
into v. Such an equation obviously conserves mass. Rotational symmetry can 
be enforced as a consequence of the assumption k{Rv* — > Rv) — k{v* — > u), for 
all V, v* , and for all rotation matrices R. Detailed balance would mean that the 
integrand vanishes, whenever / is a local equilibrium distribution. We shall only 
require that L 9JT = 0. In the right hand side / can be replaced by (1 — 11) /. It has 
been shown in T that in this case an H-theorem holds: 

^ m. 4 iRd V 2)T m* J \m m*J 

where /* denotes f{v*). A sufficient condition for microscopic coercivity is 

Assumption (Hl.l) ^^^ + tt^ > 2A„. > 0. 

Note that, because oi LdJl — 0, the collision frequency i^(v) = J^^ k{v -^ v*) dv* 
can be written as 

Thus, (Hl.l) implies z^(v) > Am. 
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For proving the boundedness of A L, note that g — ALf can be written as 
g — u F with 



ue 



V^ • (e-^V:,^) 



x-jLf 



Muhiphcation by u and integration gives 



/ u^e dx + |Vj;Up e dx — VxU-j\_fdx. 

For the gain term L+/ = L^ k{v* — ?■ v)f* dv*, by applying the Cauchy-Schwarz 
inequahty twice, we get 



\J 



L+/I 



< 



< 



\r 



m* 



\f* 



m* 



■ dv* 



dv* / k{v* -^ vym* dv* 



m* 



dv 



\v\ fc(w* — > w) '7^ dv* dv 



An analogous estimate for the loss term holds, so that we finally have 

y.u-j^fdx <c(fjVxufe-^dx] 11/11 
under the assumption 
Assumption (H4.2) // {\v\^ + \v*\^) k {v* -^ vf -:zr dv* dv < oo . 

1/2 



As a consequence, ||AL/|| = (/jj^lupe ^ dx) 



|2„-y^^^-/' < (^lljll aj^j AL js bounded. 
Combined with (H4.1), (H4.2) shows that (H4) holds. 

Theorem 11. Let L be given by ([32]). // (HO.l), (Hl.l), (H2.1), (H4.1) and 
(H4.2) hold, then solutions of (jl.ip with initial data in L^{dfJ,) decay exponentially 
to the global equilibrium given by (|3.1[) . 

4. Linearized BGK operators 

4.1. Motivation: nonlinear models. Our motivation in this section comes from 
nonlinear BGK models with collision operators of the form 

1 



The operator is determined by the energy profile ^(E) > which is assumed to be 
monotone decaying on 7~^(0, oo). The (strictly increasing) function /J(p) is defined 
implicitly by the requirement of local mass conservation, i.e. 

7 {t:^^ -'Pip)] dv = p. 

Global equilibria of the nonlinear equation dtf+T f = Q{f) are given by fooix, v) = 
'y{E{x,v) — fJ-oo), where the constant ^cxd is determined by the total mass, and the 
macroscopic equilibrium density poo by p^o — V{x) = 'p{poo{x)). In this section, 
we shall investigate the linearized stability of these equilibria, leading to the linear 
equation (jl.ip with the linearized collision operator L = 11 — 1 with 

F{x,v) = -i{E{x,v)- p^), Pf(x) = ^Tf TTY- 

M [Poo[X)) 
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Note that in the Maxwelhan case ^{E) — e^^, Ji{p) = log(p(2 tt)"''/^), the operator 
Q is hnear, and therefore equal to L. This case has aheady been investigated in 
Section 13.11 

The macroscopic hmit 

dtp^y^- (V^i^(p)+pV^y), 

of the nonhnear equation is a drift-diffusion equation with nonhnear diffusivity 
cr(p) — v' {p) — p'p!{p) (see |10| for a justification). Macroscopic limit and lineariza- 
tion commute in the sense that the linearization 

dtp = ^x ■ (Va;(cr(/9oo)p) + P^xV) 

of the macroscopic equation is the macroscopic limit of the linearized kinetic equa- 
tion. 

In the following section we consider a family of equilibrium energy distribu- 
tions 7, giving rise to nonlinear diffusions of fast diffusion type. As in Section [3. 11 
boundedness of L and microscopic coercivity are straightforward. Since j\_f = — j/, 
A L is bounded, and L A = is easy to check. It remains to check the macroscopic 
coercivity condition and the boundedness of AT (1 — 11), corresponding to (H2) 
and (H4) respectively. 

4.2. Fast difTusion. The choice -i{E) = £;-'i/2-i/(i-m) and m < 1 leads to 

, , _ r p™ for m 7^ , 
[ log p for m = , 

with a constant c > depending on m and d. For ?7i < 1, we also compute the 
moments 

PF - Co(y - Moo)-'"'/^'""^ , TOF = Ci(y - poo)-^'^^-^^ , 

MF^C2iv-p^y-'/^'-'^\ 

where the positive constants cq , Ci , C2 depend on m and d. For the external potential 
we shall, for notational convenience, only consider the choice 

(4.1) V{x)-p^ = {l + \x\y, (3>0. 

However, all our results are easily extendable to potentials whose asymptotic be- 
haviour as \x\ — ?► oo is given by (|4.1|) . With these choices, 

d(l — m) 
Assumption (HO. 2) /3 > —7- '-. 

2 (2 — m) 

is necessary and sufficient for pp e L^{dx). 

Macroscopic coercivity is related to Hardy-Poincare inequalities. In [11131 13|, for 
any d > 3, a 7^ a* := —{d — 2)/2 a positive constant Ca.d is given explicitly, such 
that 



(4.2) / \Vxumi + \x\^rdx>Cc.^d [ m'(1 + 



^r-Ux. 



for allu e H^{{l+\x\'^) dx), under the additional condition Jj^^ u (l + |.Tp)"^^ dx = 
if a < a*, in which case the measure (1 -I- |a;p)"~^ dx is bounded. The Hardy- 
Poincare inequality is equivalent to macroscopic coercivity (H2) for /3 = 1. A small 
generalization is even more useful for our purposes: 
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Corollary 12. Let d > 3, ai > a2 + 1, and ai ^ a* := 1 — d/2 if ai = a2 + 1. 
Let w be a function such that < w{x) < c (1 + |a;p)"^ for any x G W^ , for some 
c > 0. Then there exists a positive constant ICai.a2.d such that 

f |V,up(l + |xn"ida;>/Ca,,„,,d / u^wdx 

for any u such that J^^^ uw dx — if J^^^ w dx < oo. 

Proof. The assumptions on ai and a2 allow to choose a ^ a* with a2 + i < a < cti- 
Then Theorem 1 in [2 implies 

\V,uf{l + \x\^)"'dx> f |V,u|2(l + |x|2)"dx 

>C,d/ \u^\Hl + \x\^r~'dx>Ca,,d f \u^\^wdx, 



where u = u + w^ with u = for a > a* , and u = Const and J„a u^(l + 
\x\'^)°'~'^dx — for a < a* . This completes the proof for a > a* . Otherwise, 
/jjji wdx < oo holds, and the right hand side can be estimated as follows: 



|u — u| w dx > inf / \u ~ fi\ w dx ^ u wdx, 

MeRjRd jRd 

using the side condition J^^^ uw dx = 0. D 

To get examples, where macroscopic coercivity holds in the fast diffusion case, 
we apply Corollary [T2l with ai = — /3/(l — rn) and a2 = — /3 (2 — w)/(l — m). Then 

Assumption (H2.2) d > 3, /3 > 1, and m ^ {d - A)/(d - 2) 

implies (H2). Note that the last condition is needed only ioT /3 — 1. It will however 
be useful in the following. For proving the boundedncss of AT (1 — H), a modified 
version of the framework of Section [2] can be used. Redefining 

'^ "'"''- < cW holds with the notation of Section [2] The result of Lemma [T] 



(4.3) \\WV,u\\1>k' 



W^u 



is a direct consequence of the Hardy- Poincare inequality with a = f3 — l + l3/{m—l). 
Note that we have to require that a is different from a* . The proof of Lemma[5]uses 
Assumption p.Sp . which would together with (H2.2) require /3 — I and therefore 
by (HO. 2) m > {d — 4)/{d— 2). This can be slightly improved by redoing the proof 
of Lemma [51 

Lemma 13. Let (H2.2) hold. Then there exists a constant /3q > 1, depending on 
m and d, such that for (i < j3q the operator AT (1 — 11) is hounded. 

Proof. As mentioned above, (|4.3p follows from [2] Theorem 1]. Since, by (H2.2), 
l/(m — 1) ^ a*, also /3 — 1 + /3/(m — 1) 7^ a*, if (i is close enough to 1. According 
to [1], the explicit expression of the constant Ca.d in (|4.2I) is a positive, continuous 
function of a for a < a* and a> a* and n' — k'(/3, to, d) > in ()4.3p can be chosen 
to be continuous with respect to /3 at /3 = 1. 

As in the proof of Lemma HI we derive the inequality (|2.10p and use 

wl |V.(W^')| < 2 (/3 - 1) wo wi VF^ 
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to estimate the last term by 2 (/3 - 1) ||m VF^Ho ||P^ V:eu||i. With the help of (|43| . 
(|2J0l) implies 

\\W V Mil < ^\\WV^u\\i{\\uf Wo + 2 {P-l)\\WVM\i'^ 

By the continuity of k' , there exists (3o{'ni,d) > 1 such that 2 (/3 — 1)/k' < 1 for 
1 < /3 < I3q. or such a /3 the result of Lemma [8] follows. This allows to carry out the 
proof of Proposition [SJ since Assumption (|2.6p . which is used there, is satisfied. D 

As a consequence of this result, we formulate 
Assumption (H4.3) /3 < /3o with /3o("i,(i) from Lemma [T3l 

Theorem 14. With the above notations, let L = II — I, and assume that (HO. 2), 
(H2.2), and (H4.3) hold. Then solutions of (jl.ll) with initial data in L^{d^) decay 
exponentially to the global equilibrium given by 

/-, X -d/2-l/(l-m)-l 

F{x,v)=[^-\vf + V{x)j , V{x) = {l + \x\y, 

This result is, to our knowledge, the first hypocoercivity result for kinetic equa- 
tion whose Gibbs state does not separate position and velocity variables. 
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